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Abstract 

In a new functional integral approach proposed for the model, we find 
the regime with a deformed integration measure in which the standard in- 
tegral is replaced with the Jackson integral. We indicate the relation to a 
p-adic functional integral. For the magnetic and electronic subsystems in 
the effective functional that results from the operator formulation of the 
Hubbard model, we find the two-parametric quantum derivative result- 
ing in the appearance of the quantum SUrq{2) group. We establish the 
relation to the one-parametric quantum derivative and to the standard 
derivative. 

1 Introduction 

The Hubbard model |T] remains the main object on which to probe new ap- 
proaches to strongly correlated systems. It is traditionally assumed that ef- 
fects of strong Coulomb repulsion play an important role in understanding 
mechanisms of high-temperature superconductivity, the physics of a metal- 
dielectric transition of the Mott-Hubbard type, and the related magnetic states 
[2]. The experimental discovery of spin liquid in an organic crystal {BEDT — 
TTF)2Cu2{C N)^ [3] made the situation even more complicated, establishing 
the necessity of taking spin liquid into account in the phase diagram of the 
model under study. It recently became clear that we need new approaches for 
help, in particular, in understanding properties manifested by materials and 
separate molecules used in nanotechnology devices. Materials with a one- or 
two-dimensional conductivity band, for instance, polymers, are widely used in 
these devices. For the TA-PPE polymer (the abbreviation of poly(p-phenylene 
ethynylene), jumps of two types, step and sawtooth, were observed [7 in the 
dependence of conductivity on the applied voltage. 

Jumps of the conductivity and magnetization were observed in other systems 
and materials [H 13 16] . We proposed ||8j describing such a behavior of macro- 
scopic quantities using functions of the argument taking values in the p-adic 
number field. The argument of this function determines the value of the volt- 
age, while the numerical simulation has demonstrated that the function itself. 
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Figure 1: Fig. 3. The dependence of the conductivity on the voltage for the 
polymer TA - PPE: in the insets, we present magnified jumps in the ranges (a) 
from -4.0 V to -5.0 V and (b) from 4.5 V to 5.0 V. 



which depends on a one parameter, results in the two experimentally observed 
dependence types: step and sawtooth. Increasing the accuracy when calculating 
this function of the p-adic argument (we took the values p = 2, 3 for the p-adic 
numbers) revealed the appearance of a nested systems of jumps in the conduc- 
tivity behavior. The physical pattern of such a behavior in strongly correlated 
systems was described in [3] on the qualitative level. It was mentioned that in 
the regime of strong electron repulsion in one-dimensional systems, the electron 
system splits into fluctuating clusters containing several electrons each. As the 
temperature and voltage increase, the viscous motion of these clusters taken as 
a whole results in processes of splitting them into smaller clusters. We thus ob- 
tain a hierarchical embedding of jumps of smaller amplitude into the structure 
of jumps with larger amplitude. 

This paper is devoted to justifying that a regime in which the studied macro- 
scopic quantities are described by functions of a p-adic argument can arise in 
strongly correlated multielectron systems with strong Coulomb repulsion. We 
here justify the following scenario: p-adic numbers appear through a deforma- 
tion of the integration measure depending on the deformation parameter q. The 
Jackson measure used when calculating macroscopic quantities appears in the 
functional integral in the strong correlation regime. In the case q = where 
p is a prime, the Jackson integral becomes the p-adic integral |18j . We can 
therefore use functions of a p-adic argument to describe jumps in the Hubbard 
model. In particular, we here continue the investigation of the new formulation 
of the functional integral that we proposed in |16} [TT] . We note that this ap- 
proach differs substantially from functional formulations of the Heisenberg and 
Hubbard models common in the foreign literature. Therefore, in what follows, 
we present a comparison with different versions of functional representations for 
creation-annihilation operators proposed in various papers. We show that the 
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obtained expressions for supercoherent states result in more complex composite 
expressions for the operator symbols than those provided by the slave-boson and 
slave-roton approaches |14l fT5] . The expressions for the creation-annihilation 
operators obtained here with various group reduction schemes in the Hubbard 
model first proposed in |16| . result in all currently known formulations of the 
functional integral for systems with strong Coulomb interaction. Moreover, as 
shown in |T7j, the proposed formulation allows studying various cohomologies 
of groups and supergroups and thus provides a controlled concretization and 
expansion of dynamical groups and symmetries spontaneously appearing when 
strengthening the interaction in multielectron systems. Our proposed represen- 
tation for the effective functional and creation-annihilation operators is essen- 
tially nonlinear. This nonlinearity allows segregating those terms that provide 
quantum derivatives determining generators of the quantum algebra in expres- 
sions belonging to the universal enveloping algebra. In |10| . we introduced the 
approximation in which radius vectors of vector fields were independent of the 
dynamical field coordinates. This assumption is equivalent to the approximation 
in which the values of the fields that correspond to the Casimir operators and 
to the invariants of the classical groups SU{2) and SU{1, 1) in our formulation 
are coordinate independent. Below, we demonstrate that these combinations 
determine two parameters of deformation of the quantum Lorentz group, which 
is the group of four-dimensional rotations. We treat the transition to the one- 
parameter quantum derivative in detail and show that the magnetic and electric 
subsystems of the Hubbard model are described by deformed versions of nonlin- 
ear sigma models. We further calculate the part of the functional integral for the 
Hubbard model that results from expressions containing this quantum deriva- 
tive. We realize our consideration on model contributions naturally present in 
the Hubbard model. We consider the limiting cases with respect to the defor- 
mation parameter and calculate the functional integral in these limit regimes. 
We calculate the effective functional following from the kinetic energy of the 
Hubbard model using supercoherent states containing no more than the first 
powers of fermionic fields. We fix the dynamical fields such that the residual 
expressions indicate the difference between our approach and other approaches 
in the literature. The main distinct feature of our representation is its non- 
linearity, which allows introducing a two-parameter quantum derivative into 
the problem. When passing to a one-parameter quantum derivative, we obtain 
a deformed nonlinear sigma model in the auxiliary dimension determined by 
the scale factor. Further, by analogy with the diagram technique for atomic 
X-operators [21 [IB], we consider the integral of the scale factor linear in the 
quantum number. This integral determines the integration measure. Directly 
calculating the corresponding series, we show that in the limit case as the de- 
formation parameter tends to zero, this series converges to the Jackson integral. 
Obtaining deformations of the integration measure when varying the deforma- 
tion parameter is our main result. In the limit case of small q, the standard 
integral becomes the Jackson integral. 



2 The effective functional of the Hubbard model 

We briefly describe the general scheme for constructing the effective functional 
for the Hubbard model proposed in [TU] and recently further developed, for 
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example, in [TT]. We begin with the Hubbard model written in terms of the 
standard creation-annihilation operators, 



H 



mE^ 



(1) 



where Q;j^j,Q;a-j are the creation and annihilation operators, n^^i is the elec- 
tron density operator, W, U, jj, are the respective width of the conductivity band, 
the one-site repulsion of two electrons, and the chemical potential, and a deter- 
mines the spin value. We sum over the indices i and j labeling the lattice sites. 
For instance, for a one-dimensional lattice, we have ~N/2 < i < N/2, where N 
is the total number of atoms. For electrons, we sum over the indices a — ±1/2. 
In the weak correlation regime with a small on-site repulsion, we take the zeroth 
approximation for the kinetic energy and treat the repulsion perturbatively. In 
the "atomic" approach with strong repulsion, we take a one-site contribution of 
the form 

as the zeroth approximation. In what follows, we omit the index i in order to 
concentrate on the internal structure of the appearing fiber bundle. This term 
is diagonal in the "atomic" basis and has the eigenfunctions and eigenvalues e: 

|0 ^; 1+ a+|0 |- a^|0 ^; |2 ^= a^a^|0 ^, 



"l^ I 

0,e+ 



Calculating the matrix elements for a'^^a^ in this basis, we obtain the fol- 
lowing matrix representation. For instance, is: 
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(2) 



The operators a'^ therefore contain two nonzero matrix elements. We introduce 
the operators X^'^ , r,s = 0,+,—,2 containing only one nonzero matrix element. 
We then obtain the expansion of the creation and annihilation operators in 
terms of the Hubbard operators [T]: 



a+ = X+" - X^-,a+ = + X^+,a^. = X"+ - X'^, = X°- + X+\ 

Because we have a four-dimensional basis and the creation and annihilation 
operators are expressed as 4x4 matrices, we have a basis comprising 16 opera- 
tors in the general case. From this set we remove the unity operators and the 
operators of the form 



75 = 
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We use this notation for this operator because it coincides with the 75 op- 
erator written in the chiral basis in quantum field theory (QFT). It also equals 



75 = - x^^f - {X++ - x—f. 

We can express the one-site Hubbard repulsion in terms this operator. Be- 
cause we can express it in terms of other operators, we do not take it into 
account in what follows. The remaining operators on a site in the given basis 
can then be separated into the fermionic operators of the form 

XO-, X-\X+\X-\ ^2+, 

and the bosonic operators of the form 

{X+-,X-+,X++ - X— , ^20, - ^22). 

We omit the lattice index of these operators to indicate the possibility of ob- 
taining generators of the global dynamical algebra. We can endow the Hubbard 
operators with the lattice index by taking N their copies, where N is the total 
number of lattice sites, and constituting the direct product from these copies. 
A contemporary review of the atomic approach was presented in In terms 
of these operators, the Hubbard model becomes 

H^uY^xr ^wY, xrx^ (3) 

In what follows, we use the functional formulation for multielectron systems 
proposed in [TU]. The evolution operator between the initial and final states 
in this formulation is given by the following functional integral with the action 
expressed in terms of the effective functional, which is to be calculated using 
operator expression ([s]) using the supercoherent states (we note that the index 
labeling the lattice in this formulation is lacking for the set of operatorsX""* , 
which therefore constitute the set of generators of the dynamical superalgebra) : 



AG(tf)> = \Gf> 

<Gf\e-'"^*f-'^^G,>= I?(G,G*)e-'^[^'<^ 1; (4) 

J\G{ti)> = \G,> 



The action here has the form 



<G{r,t)\i§-^-H\G{r,t> 



sican^l ,tj^,'r --,^(^-- (5) 

and the integration measure is 

D(G,G*)= Y[ Y[dG{r,t)*dG{r,t) <G{r,t)\G{r,t)> 

ti<t<tf r 

The integrand in QFT formula ([s]) is called the Lagrangian. We avoid using 
this term wherever possible because in contrast to QFT, we begin with an al- 
ready secondarily quantized representation and subsequently pass to an effective 
model containing composite and nonlinear fields. Such a passage is nontrivial, 
and to use the term "Lagrangian," it is desirable to know the Lagrangian equa- 
tions of motion at least. Unfortunately, in our approach, for each representation 
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of the chosen supercoherent state, we must separately investigate the search for 
the variables and their conjugate momenta. We therefore use the term "effective 
functional" (although it is not very precise) instead of the term "Lagrangian." 
We use a function equal to the supercoherent state of the form 



I G >= exp 













Xi 




h+ 





X2 


h- 


-K 







~X3 


X4 





to describe local properties of the strongly correlated system. The exponent 
in this expression contains dynamical fields depending on the coordinates and 
time. We determine the electric field by the three-dimensional vector, 



E = {E+{x, y, z,t),E {x, y, zt), E'{x, y, z, t)). 
magnetic fields has three components depending on the space-time ordinates, 

h = {h'^{x,y,z,t),h^{x,y,z,t),h^x,y,z,t)). 

The fermionic fields are described by odd Grassmann- valued functions of the 
coordinates and time, 

xUx, y, z, t),Xk{x, y, z, t), k^l,2, 3, 4. 

Expression ^ is the solution of the following problem: in the atomic basis, the 
site state is described by a wave function of the form 

I G {i) >= a, I lO) + h I i+) + c, I 1-) + d, I z2), (7) 

where ai,bi, Ci, di - are the coefficients of the atomic basis functions depending 
on the index i. We have the normalizing condition for the wave function 

< G I G >2= +b^ + c^+dj = l 

The drawbacks of using functions of this type are well known. A function 
of the typeci is symmetric with respect to its coordinates (bosonic type) , while 
the functions bi , Ci are skew-symmetric under coordinate interchanges (fermionic 
type). The general form of the functions | G > is expression where the 
exponent has the form J2mn ■^™"'^'mn,3^^^<i the functions are Xut^tE. Su- 
percoherent state (|6| determines the supergroup representation with the set of 
supergenerators X"™. We note that the expressions for the above supercoher- 
ent states were calculated exactly in [20]. As noted in [17], using computers 
and computer programs for working with supermatrices and superpoles, we can 
exactly calculate the functional itself. It is clear from [20j that the main distinc- 
tion between our approach and the approaches in other papers is the nonlinear 
form of the supercoherent state expression. Below, we show that all the rep- 
resentations proposed in the literature for symbols of the creation-annihilation 
operators follow from our representation when the fields are fixed. We also show 
that it is the essentially nonlinear nature of the expressions in the supercoherent 
state, which are functions of electric and magnetic fields, that distinguishes our 
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representation from the previously proposed functional integral representations 
for the Hubbard model. These nonlinear expressions contain quantum deriva- 
tives and, as a result, quantum groups. We note that after the super- Yangian 
structure in the one-dimensional Hubbard model was discovered [T^ , there were 
numerous attempts to find quantum groups in the Hubbard model, but they 
were unsuccessful. We use an absolutely different approach in this paper, and 
the discovery of quantum structures in the Hubbard model is therefore explain- 
able. The expression for the effective functional S governs the transition from 
the operator formulation in terms of the creation-annihilation operators and the 
Hubbard operators to the field theory formulation in terms of dynamical fields of 
the Bose and Fermi types. As shown in [17], the main point of this transition is 
the passage to the problem formulation in terms of the superbundle determined 
by a representation of the local supergroup of four-dimensional space rotations. 
In what follows, we need only the symbols of the creation and annihilation op- 
erators determined by the expressions < G \ \ G>,<G\aa\G>. The 
general expressions for supercoherent states were calculated in [SU]. Here, we 
take expressions for these states in the zeroth- and first-order approximations 
with respect to the Grassmann fields. We write the expressions for the coherent 
states and their conjugates in these approximations: 

/ ch{E)+E/-^ \ 
I Q a+xi + atx2 

aiXi+a2X2 



<G\^icMiE)-E^^,iatrxl + ia+rx*2A^irx*i + ia^rx*2,-E+^) 
The action of the operators on this state leads to the expressions 



a+ \G>= 



/ ^^ 

chiE)+E,'-^^ 



,a+\G>^ 



G >-- 



\ -«! Xi -0-2X2 J 
/ atxi + atx2 \ 




/ 


cME) + E,'-^^ 
V a+xi + ajx2 J 



-E 



V 



_ sh(E} 
E 





/ a^Xi+ 0-2 X2 \ 



,ai\G>= 



J 



E 



- sh(E) 



\ 



E 



J 



The symbols for the creation and annihilation operators calculated on these 
states and written in the spinor form are 



<G\a+\G> 
<G\a+ \ G> 

<G\a^\G> 
<G\a^\G> 



= E^20mX + EllAmX*, 

= E22amx + E2iAmx*- 



(8) 



Here, the 2x2 matrices E and hi are 
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Ei2 = £'12 ( ^ ^ j , E21 = E21 f 0^ 1 



Ell E12 \ _ f ch{E) + E,_sh{e)/E E+sh{E)/E 
E21 E22 ) ~ \ E-sh{E)/E ch{E) - E^sh{E)/E 



dm = h{Hhlf2))' + h{{Hf2))' + E,)h2. 

{In{f2)y + K h+ 

h- {ln{f2)y - K 



h. = 



1 

1 



We recall the fields over which we integrate in the functional integral: the 
vectors E and h describing fluctuations of the respective electric and mag- 
netic degrees of freedom and the Grassmann fields Xi- The combinations E ~ 
y^E'^ + E+E^ , h — \/h1 + h+h^ , which are invariants of the two SU{2) groups, 
also enter the supercoherent state. We here set them to be coordinate- and time- 
independent constants. These constants determine two deformation parameters 
in the quantum group problem. We now consider the coefficients /2 and /a. Ex- 
pressions for these quantities are nonlinear in the above invariants. The method 
for calculating them exactly was developed in [20 . In the nonlinear represen- 
tation of supercoherent states, we have four such coeficients:/, /2, /a, /4. We 
describe the trick used to calculate these coefficients in [17]. We multiply the 
invariants E and h in the formulas by the scale factor x and subsequent show 
that we can obtain /, /2,/3,/4 from / by taking derivatives with respect to 
X. At the end of the calculations we set the parameter x equal to unity. In 
the formulas presented above, such differentiations are marked by the prime. 
Expression ([8| determines the representation for the creation-annihilation op- 
erators in the functional integral in terms of the three components of the vector 
h. The representation depending on the vector h was proposed in [14 . We note 
that in the expressions there and also here, the local group is SU{2)* SU{2). We 
consider particular cases of the above composite expressions for the creation- 
annihilation operators: 1) If {E^ , E^ , E^) — 0, {h'^ , h~ ,hz) — ,then: < G \ 
01^ I G >= X* , < G \ \ G >~ Xa we obtain the holomorphic representa- 
tion for the creation-annihilation operators. Substituting them in Q, we obtain 
the standard representation for ferminonic systems in terms of odd Grassmann- 
valued fields [13 2) If /i", /i^) = 0, {E+^E'^E^) ^ then expression ^ 

generates the following formulas for the functional symbols of the creation and 
annihilation operators: 



< G I a+ I G > 

< G I a+ I G > 



E'i2hBx + EiihB*ax*, (9) 



<G\a^\G> \ _ 
<G\a^\G> 



= E22hBx + E2ihB*ax* (10) 



where B — {ln{f-i))' + E^- Setting the expression for //a equal to unity, we 
reduce these expressions to the slave-roton representation in |15| . 



8 



3) If {h^,h~ , hz) 7^ 0, {E~^ , E~ , Ez) = , then expression (8) generates the 
following formulas for the functional symbols of the creation and annihilation 
operators ^ 15]: 



<G\a+ \ G> 
<G\a+ \ G> 

<G\a^\G> 
<G\a^\G> 



= E,iAmX*, (11) 



= E22amX (12) 



0-m — h^i E^nd in the logarithmic function is f^, in this formula, we have 
A^a — cro*m- We have thus shown that each representation for the creation- 
annihilation operators in the Hubbard model known to us is contained in some 
particular case of expression (([8| . We note that the main feature distinguishing 
our representation from those in the literature is the presence of the coefficients 
/2 , /a H /4- Expressions for these coefficients were given in [20]. Below, we 
clarify the consequences of the presence of these coefficients. 



3 Quantum derivatives in the Hubbard model 

We see which functional results from the operator expression for the magnetic 
subsystem of the Hubbard model. Substituting expressions for the creation- 
annihilation operators in the expression for the kinetic energy, we obtain the 
functional 

We note that for //4 = 1, our representation coincides with the slave-roton 
representation used in many papers. Our aim is to study those terms in the 
functionals that are absent in the cited papers. We take the unit matrix as the 
matrix field E. which indicates the independence on the space coordinates. We 
mentioned in the introduction that we want to obtain the expression for the 
coherent state for such E. Because x is a fermionic field, when passing to the 
momentum representation, we obtain a half-filled Fermi band over which we 
integrate to obtain the total number of fermions, i.e., a constant. We substitute 
this constant for the quadratic combination of fields xtXj- the result, we 
obtain the expression for the effective functional at coinciding arguments [i = j) 

H={Uf (13) 

where the coefficient is [17 : fi = {Esh{E) - hsh{h))/{E'^ - h'^). This is 
precisely the factor that distinguishes our functional from the one in |14l [15] . 
At first glance, it is just a nonlinear coefficient. But...! We introduce the two- 
parameter derivative acting on a function f{x) by the formula 

Drgfix) = ifirx) - fiqx))/{{r - q)x) 

It is known [19 , that the three operators {Drq, x, x-^) constitute the two- 
parameter quantum algebra SUrq{2) with the following commutation relations 
for the generators {S'^ , S~ , Sz) = {Drq,x,x-^) : 
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[Sz , S^] — 5^ , [5*2 , S" ] — S , S ^ S — [215*2] 
We have the definition of the two-parameter quantum number |19j : 

[X]rq = 

Using the expression for Drq we can rewrite the nonhnear coefficient in 
the form 



DE2M^{x){y/3:shy/^/x) = {E^/^sh{^/^E) - hy/^sh{y/3:h))/{x{E'^ - h'^)). 

We introduced the separate parameter x into the expression for /4. The 
representation of this quantity in [20^ imphes that x determines the scale factor 
and is the dilaton field. We now show this. The complete set of fields in 
expression ([6] for the supercoherent state and the functional integral integration 
measure is determined by the superalgebra whose matrix generators enter (|6] 
It does not contain the dilation generator, whose matrix representation is 75 
in the chiral representation related to the third diagonal term in ([6j To take 
this term into account, we must calculate matrix (6) exactly with the full set 
of generators. Unfortunately, this problem is much harder than the problem 
solved in[2n]- We introduce the scaling factor for the local fields E and h such 
that E"— > /3E, h"— > /3h and pass to the primed fields. This factor can be 
a local x-dependent field P{x), ) or a global x-independent field. The field j3 is 
complementary to the fields E and h and must enter the nonlinear expressions 
for the supercoherent state. We can introduce a functional integration over 
this field, and this field then enters the integration measure. Because fields of 
this type (of the dilation field in QFT) had never been previously discussed 
in the Hubbard model framework, we below concentrate on discussing possible 
qualitative effects resulting from taking this field into account. In this paper, we 
set the scaling factor independently of the coordinate x. We have already shown 
that the quantum derivatives act on precisely this variable. We note that P{x), 



is replaced with x in formula ( 13 1 in order to introduce a unified notation for 
the entity of the extended space-time coordinates. Taking the scaling factor into 
account means passing to the {d + /)-dimensional space-time in the case of the 
d-dimensional initial space-time. Substituting the derivative Drq in expression 



( 13 1 , we obtain the functional of the form 

H = Y^{Drq{x)f{x)f (14) 

X 

where f{x) is given by expression in the previous formula on which the 
derivative acts. The difference from the nonlinear sigma model is due to specific 
properties of the coordinate x, which defines the dilaton field. We consider 
the limit as rr ^ (7 or 1/h^. It is clear that Drq Dq in this limit 

and becomes the one-parameter quantum derivative acting on functions by the 
formula 

Dqfix) = ifiqx) - f{q-'x))/{{q - q-^)x). 

The commutation relations for the operators ({Dq, x, x-^) constitute the quan- 
tum group SUq{2) [22 with the commutation relations 
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[Sz, S~^] — , [Sz, S ] — S , [S^ , S ] — 7/1 / \\ • 

sh[ln{q)) 

Letting g — > 1, we reduce the action of the quantum derivative Dq on arbitrary 
functions f{x) to the action of the standard derivative [21], i.e., Dq — )■ . 
Replacing Dq — )■ we obtain the standard nonhnear sigma model ■ 
We note that / is a solution of the normalizing constraint condition on the 
coherent state, which customarily defines a vector field of unit length. Such a 
model describes a magnetic subsystem in a functional integral approach. On 
the operator level, it is the Heisenberg magnet in the Hubbard model. This 
magnet model is defined on the dilaton field and is one-dimensional. Returning 
to the general effective functional in the Hubbard model, we note that we have 
found a quantum symmetry and the way to introduce it into strongly correlated 
systems. It is especially important that subsequent contractions of quantum 
derivatives coincide with three levels of the tower of symmetries discovered in 
the Hubbard model in jTB]. Formula ( 14 1, implies that in multielectron systems, 
we must study "quantum" magnets whose variables are generators of quantum 
Hopf algebras, not simple Lie algebras. We very briefly indicate the differences in 
the physical behaviors of such magnets. We here discuss only the behavior of the 
magnetic momenta and discuss more serious differences later. It is known that 
the Casimir operator Si + S'^S~ of the group SU{2) take the value: S{S + 1). 
The Casimir operator of the group SUq{2) is S^S^ + (£M|^)2 ^^^-^^ --^-^ close 
to zero as q— > 1. In the physical language, this means that at the phase 
transition point at which a nonzero value of the parameter q appears, spin 
can arise smoothly with its value changing continuously from to 1/2. This 
picture of the Mott-Hubbard-type metal-dielectric transition seems more natural 
compared with the phase transition scenario in which the spin instantly jumps 
from (metal) to 1/2 (insulator of the Mott-Hubbard type). We now pass to 
studying a qualitative a behavior of a functional of type ( 14 1 using the functional 
integral. 



4 The Jackson integral 

The above analysis of nonlinear expressions in the supercoherent state demon- 
strates that we can write complicated expressions in terms the action of the 
quantum derivative containing one or two deformation parameters. These deriva- 
tives act on functions of the variable, which was introduced in |20| as a purely 
technical tool ensuring a uniform form of writing expressions for /2, /s, /4- The 
same variable was again used purely technically in |17| to calculate and write 
nonlinear expressions. In the preceding section, we showed that by taking the 
scaling factor into account, we obtain a model resembling a nonlinear sigma 
model but with quantum derivatives instead of the standard ones. In the 
limit where the one-parameter quantum derivative transforms into the standard 
derivative, we obtain the standard nonlinear sigma model. We thus qualita- 
tively obtained a fact that is well known in the Hubbard model: the magnetic 
properties of the Mott-Hubbard dielectric at large values of the Hubbard re- 
pulsion U are described by the Heisenberg model with the exchange integral 
equal to W'^ /U [5]. Below, we use the diagram technique for spin and Hub- 
bard operators and also our diagram technique for constructing the "tower of 
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symmetries" [[H]. When constructing the diagram technique for the Heisen- 
berg model, we found the following trick useful. To the model action, which 
is quadratic in the spin operators, we add a term that is linear in the spin op- 
erators and proportional to Sz, the so-called Zeeman contribution. We then 
study the quadratic terms perturbatively. Summing the diagrams responsible 
for the local ordering results in a nonzero contribution to the magnetic field 
and ensures the possibility of calculating means over the ground state when 
setting the previously introduced field to zero. Our scheme and expression(|14|. 



imply that we encounter exactly the same situation. We cannot take ( 14 1 , as 
the leading approximation in contrast to the case where we have the standard, 
not quantum, derivative. We must add a linear derivative and study how it af- 
fects calculating the ground state means. Below, we show that by deforming the 
derivative, we do not introduce new interaction types, as could be expected, but 
rather deform the functional integration measure, i.e., we change the type of the 
functional variable itself. In addition to the standard bosonic fields described 
by complex- valued functions and fermionic fields described by odd Grassmann- 
valued functions variables, fields with a nontraditional integration measure, for 
instance, taking values in the field of p-adic numbers, appear in the functional 
integral. The base of this field is determined by the means of the bosonic electric 
and magnetic fields, which provide the deformation parameters for the quantum 
groups. The standard scheme for studying functional integral ^ begins with 
the "zeroth" approximation obtained from an effective functional quadratic in 
either bosonic or fermionic variables, which reduces to Gaussian integrals over 
complex variables parameterized by momenta. We must be able to evaluate 
func-tionals containing quantum derivatives. For subsequent use, we briefly de- 
scribe the scheme for constructing an effective functional using a holomorphic 
representation for creation and annihilation operators. To define an integral over 
trajectories completely, we must construct a functional series, which can then 
become a formal perturbation series. For this, we must know how to calculate 
integrals of the form 



x''e-'"''dx (15) 



for an arbitrary function f{x). Clearly, we can exactly evaluate Gaussian in- 
tegrals with power-law dependences, which allows standardly constructing a 
formal perturbation series. Using this, for each bosonic quantum system with 
a Lagrangian quadratic in the fields, we can set the complex numbers ip, ip* 
into correspondence with the creation-annihilation operators , /3o- , satisfying 
the standard commutation relations. Each functional integral with a bosonic 
Lagrangian containing interactions of fields can then be constructed as a formal 



series in the interaction constant. To study effective functionals of type (14 1 



we must replace the exponential in expression ( 15 1 as 



exp{bipip*) exp{bxDq) 

Because we replace the pair of operators (x, g|) with the holomorphic co- 
ordinates (if^ip*), we obtain our variant of quantization if we substitute the 
quantum derivative Dq{x). for the standard one. Prom the definition of Dq, we 
have 
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jj sh{\n{q)x 



s 



h{lTi{q))x X 



Substituting this expression in the exponent results in the expression 

/(a;)e''[^l«da; (16) 



where [x]q = ''g_g-i = '^sh^\n{q)) ■ Calculating a functional integral containing 
general quantum derivatives Drq, in a simpler case of the quantum derivative 



Dq, therefore requires calculating integrals of form (161, with arbitrary functions 
f{x) of the variable x. We have thus described the scheme arising in the Hub- 
bard model for contributions to the functional integral that contain quantum 
derivatives. Such integrals are always present in the general functional integral 



for the Hubbard model. To calculate (16 1, in the general form, we expand it in 



the Taylor series in the derivative Dq and expand the exponential [x]q in the 
series |2lj : 



c 

/. 



. n n I \Q ' 'J 



1 I 

m=On=0 







Here (x — c)q"^ — {x — c){x — qc)(x — q^c) {x — q™ ^c), and the quantum 

factorial is 

[n]q\ = [n]q[n-l]q....[2]q[l]q[0]q 



The constant c, which determines the limit scaling in the problem, arises in ( 17 1. 
We represent the action of an arbitrary power of Dq on f{x) as a series, 

n 

D'q-{x)fix) = J2 dll{q)f{q^''x) 

k——n 

A similar formula describes the action of an odd power of the quantum derivative 



on a function; the summation then ranges odd powers. Integral (17 1 can then 
be expressed as a series in /(eg™) of the form (we present only the first terms 

of the series) /(c)[/o + + j^h + + ■■■] + f{cq)[j^h + + 



+ ...] + ficq-^)[§^h + §^h + + ...]- 



Here, the quantum number [x]q = "^qZq-i tend to a; as q— > 1. We have 
the representation 
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In ^ y ~ eyi-g{b[x\q)dx 



c 

for the quantities /„. These integrals reduce to integrals of the form J cc™ exp(6[a;]g)da; 



and can be evaluated as infinite series in h (the detailed scheme for calculating 
integrals of this type and detailed expressions for the series in these integrals 
will be published elsewhere). As a result, we have /„ « 1 as g— > 0. We 
present expressions for several quantum numbers: = - 



= 1; 



= q 



q + q 
2 , i_ 



-q ^ We note that ( 17 1 reduces to ( 15 



as g — > 1. Ex- 
pression ( 17l therefore contains a standard perturbative scheme for determining 



a functional integral. We consider the opposite limit q « 1. The expressions 
for quantum numbers and their factorials in this limit are [l]g ^ 1 , [2]q ^ q~^ , 



[3], 



[4], 



[5], 



(n 1) ^Yje general term. 



Expressions for the factorials for q << 1 , are: [l]g! = 1 , [2]^ 



[3],! 



[4],! 



[5], 



and: 



.^i for the 



general term. Polynomials in the small-q regime are :p\ — —q; p\ ~ q; 
1; p| = l/g^; P5 — — 1/<Z^; Pe — 1/9^ • Substituting these expressions for q << 1 
in formula (17 1, we obtain the series 



/(a;)e^M'dx ^ /(c) + f{cq)q + f{cq')q' + f{cq^)q^ + f{cq^)q^ + f{cq')q' + . 



We find that expression (17 1 becomes the Jackson integral |2I] in the leading 
approximation in q. It is known |18| . that for q — , where p is a prime, the 
Jackson integral becomes the p-adic integral. 



5 Conclusion 

We have shown the existence of the regime described by the functional integral 
with the Jackson measure and deformation parameter q in the Hubbard model. 
At g = 1/p , in accordance with (TB], this regime is described by the p-adic 
integral. To demonstrate the features of this approach, we consider the 

following problem. As proposed in we describe the conductivity by the 

function /h(r) = J2o O'kP^'^l Ofc = (Oi 1; ■■■■,P ~ 1); ^ G The coefficients at 
are here equal to the coefficients in the p-adic number representation defined by 
the series rr — a^p^; — (0, 1, ■■■■,p — 1); /c € Z. We brieffy describe the 
necessary information concerning p-adic numbers: p is a prime; an arbitrary 
rational number r admits the representation r = p'^™ where n and m are 
coprime to p. The p-adic norm of a rational number is | r \p= p~^ , \ |o= 0. 
The field of the p-adic numbers Qp is the completion of the field of rational 
numbers Q with respect to the p-adic norm. We take p — The numerical 
modeling of this function for h = 0.5 and b = 1.5 is shown in the figures. On 
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10 20 30 40 



Figure 2: Fig. 1. A subcritical wave for b < 1. 

the ordinate, we have integers presented as series for r. The coefficients in fb{r) 
are taken from the p-adic expansion for these numbers. The result for fb{r) is 
presented on the axis OY. We set b = 0.5 in Fig. 1 and b = 1.5 in Fig. 2. In 
Fig. 3, we present the dependence of the conductivity on the voltage for the 
TA-PPE polymer. It is clear that the proposed function of the p-adic argument 
qualitatively describes two types of steps in this compound. 

We have shown that the regime described by the p-adic functional integral 
exists in the functional integral for the Hubbard model. 




Figure 3: Fig. 2. A supercritical wave for b > 1. 

Finally, we stress that the functional integral in the Hubbard model, as 
shown in this paper, can be successfully reformulated in terms of both the stan- 
dard and the Jackson integral and that the latter can be reduced to the p-adic 
integral in a particular case. We have the impression (supported by the compar- 
ison with experimental data [HEIEIIZ]) that in the strong correlation regime, 
the description of multi-electron systems is invariant under the choice of the 
number field. We can begin with either a p-adic or real number description. It 
seems plausible that this situation can be clarified if an exact analytic calcu- 



lation of coefficients in expressions (17 1 would be possible. Because our main 



regimes were q— > 1 and q— > 0, it seems plausible that we should study 
the crossover regime from the p-adic to the standard perturbation series based 



on Gaussian integration in more detail if we will be able to calculate series ( 17) 
more exactly. This regime would clarify an interesting question: which order 
parameter appears in the system upon the introduction of the p-adic descrip- 



15 



tion whose distinct feature is the appearance of fractal behavior in physical 
quantities? 

Acknowledgments. The author thanks the referee whose questions helped to 
improve the paper substantially. 
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